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In this paper we explore the transmission properties of single subwavelength apertures perforated 
in thin metallic films flanked by asymmetric configurations of periodic arrays of indentations. It is 
shown how the corrugation in the input side can be used to transmit selectively only two different 
wavelengths. Also, by tuning the geometrical parameters defining the corrugation of the output side, 
these two chosen wavelengths can emerge from the structure as two very narrow beams propagating 
at well-defined directions. This new ability of structured metals can be used as a base to build 
micron-sized wavelength de-multiplexers. 
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I. INTRODUCTION 

The new technological abilities for structuring met- 
als in the micro and nanometer length scales are trig- 
gering a renewed interest in the optical properties of 
metals^i&^ei^. a. breakthrough in the field was the 
experimental discovery- that the optical transmission 
through an array of sub- wavelength holes is, for certain 
wavelengths, boosted orders of magnitude over what a 
theory of independent holes would predict. 

In subsequent works it was shown how the transmis- 
sion through a single aperture can be boosted if it 
is appropriately flanked by surface corrugations in the 
side light is impinging or&ifliii. Corrugations on the 
exit surface practically do not modify the total trans- 
mittance but, for some resonant wavelengths, greatly 
modify the angular distribution of the transmitted light, 
leading to strong beaming effectsii*^. This was shown 
for a hole surrounded by corrugations in the so-called 
bull's eye geometry and for a single slit flanked by one- 
dimensional grooves. It has been shown that these two 
processes (enhanced transmission and beaming) work 
independentl y 12 ' 13 , so it is possible to have more light 
passing through an aperture and have it transmitted in 
a narrow beam by patterning both the input and exit 
sides. These effects are associated to the formation of 
electromagnetic resonances (leaky modes) running along 
the metal surface, and radiating as they move due to the 
coupling to homogeneous waves in vacuum. The spec- 
tral position of these EM resonances depend on the ge- 
ometrical parameters defining the corrugations (distance 
between indentations, and their widths and depths). 

These studies were done in a symmetric configuration 
of a finite number of grooves surrounding the central slit. 
In this paper we show that, to a large extend, the prop- 
erties of the groove arrays at the left and right of the 
slit behave independently. This property adds to the in- 
terpretation of the properties of the system in terms of 
running surface waves. It also allows, through the appro- 
priate choice of geometrical parameters for left and right 
groove arrays at both entrance and exit surfaces, this 



type of structures to act as wavelength de- multiplexers: 
the input side of the device can be used as a filter for two 
particular wavelengths and the corrugations on the exit 
side can be tailored in order to beam each of these two 
chosen wavelengths in different directions. 

The paper is organized as follows: in Section 2 we 
present the theoretical formalism. Section 3 presents 
the transmission results for a single slit with an array 
of groove to one side, while the transmission angular dis- 
tribution is considered in Section 4. In Section 5, we show 
how this kind of structures may be used for wavelength 
demultiplexing (WDM). 



II. THEORETICAL FORMALISM 

We are interested in analyzing the transmission prop- 
erties of the system depicted in Fig.l: a single slit flanked 
by finite arrays of one-dimensional rectangular inden- 
tations. The central slit may be flanked by arrays of 
grooves at both input and output metal surfaces and, in 
each case, the arrays may be different to the left and 
right of the slit. We refer to these four arrays as IL, 
IR, OL and OR (see Fig.l), depending on whether they 
are input-left, input-right, output-left and output-right, 
respectively. Each groove (a) is characterized by its posi- 
tion (r a ), its width (a a ) and its depth {h a ) 1 while the slit 
(labelled by the sub-index 0) is considered to be placed 
at the origin, having a width ao- Metal thickness is W . 

Let us concentrate in the case of p-polarized light for 
which interesting resonant effects have been reported. 
Electromagnetic (EM) fields associated to this scatter- 
ing process are expanded by a set of plane waves in vac- 
uum and by a modal expansion inside the indentations 
(slit and grooves). For this last case, as we are mainly 
interested in sub-wavelength apertures, we consider that 
only the fundamental propagating eigenmode is excited. 
Then, for x inside indentation a, the wavefield is a lin- 
ear combination of <j) a (x) exp(±zfcz), where k — 2tt/X 
and 4> a (x) = l/^/a^. We assume that the metal behaves 
as a perfect metal, expelling the electric field from its 
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FIG. 1: Schematic illustration of the structure under study. 
IL, IR, OL and OR label the input-left, input-right, output- 
left and output-right arrays of grooves, respectively. The ref- 
erence system considered is also plotted. 



bulk. One advantage of the perfect metal approximation 
is that as all length scales associated to the dielectric con- 
stant (skin depth and absorption length) disappear from 
the problem, all results are directly exportable to other 
frequency regimes (microwave or infrared regimes), sim- 
ply by scaling appropriately the geometrical parameters 
defining the structure. With respect to the applicability 
of this approach in the optical regime, we have demon- 
strated in previous theoretical worksii that this is a rea- 
sonable approximation when analyzing optical properties 
of nanostructured metals like silver or gold in simpler 
systems like reflection gratings or ID arrays of subwave- 
lcngth slits. Moreover, within this framework, we have 



been able to reproduce 
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m semi-quantitative terms 



the reported phenomena of enhanced transmission and 
beaming appearing in the optical regime for finite struc- 
tures. 

In order to extract the expansion coefficients in vacuum 
(related to the reflection and transmission functions) and 
indentations regions, we match the parallel components 
of the EM fields (E x and H y ) at the interfaces present in 
the system. By projecting the matching equation for E x 
at z — — W and z = onto plane waves, we can express 
the reflection and transmission coefficients in terms of the 
wavefield expansion coefficients inside the indentations. 
H y must be continuous only at the indentation open- 
ings and also (as we are assuming perfect metal bound- 
ary conditions) E x = at the bottom of the grooves. 
Therefore, all matching equations can be expressed as a 
function of the set of expansion coefficients inside the 
indentations: [E^,E%}(see Fig.l). The sub-set (E^) 
give the x-component of the electric field at z = —W + 
through the equation: E x (x, z = — W + ) = J2 a Ei<t>a(x), 
where a runs from — N[ to (N[ and are the 

number of grooves to the left and right of the slit in 
the input side, respectively). {Eg) is related to the 
x-component of the electric field at the exit interface 
z = 0~:E x (x,z = 0~) = ^2pEp(j)p{x), (3 running from 
-Ng to +N° (Ng and Ng are the number of grooves 
to the left and right of the slit in the output side, respec- 
tively). In both expansions, sub- index stems for the 



central slit. After some elementary algebra, it is possible 
to find the set of (N[ + + Ng + Ng + 2) equations 
that govern the behavior of [E^E^]: 



[G aa — £q] E^ + GayE^ — G V E° 8 a Q — I a (1) 

\G m - e ] E° + J2 - G V E^ = 



The first set of equations accounts for the expansion 
coefficients at the input interface and the sub-indexes a 
and 7 run from — N[ to N^. Due to the normalization 
we use for the transmittance (see below), the illumina- 
tion term I a is equal to 2 for all a. The second set of 
equations is homogeneous (the output interface is not di- 
rectly illuminated), and /3 and r] ranges from -Ng to 
Ng. Note that the only connection between the two in- 
terfaces is via the central slit and appears through the 
term GyE^SaQ in the first set and by G v EqS/3q in the 
second one. In RefsiiSii^ we explain in detail the origin 
and values of the different terms appearing in Eq.(l). For 
instance, G a pEp is the contribution to the E x field at in- 
dentation a coming from light scattered at indentation 
(3. In fact, G a p is the projection onto wavefields 4> a and 
4>I3 {G a /3 —< <j) a I G I (j)p >) of the Green's function: 
G(r,r') = fH^\k \ r - r' |), being the Hankel 

function of the first kind. G aa =< <f> a \ G \ <j) a > is 
a self-interaction that takes into account the coupling 
of the eigenmode at indentation a with the radiative 
modes in vacuum regions. This quantity depends only 
on the ratio between indentation width and wavelength 
of light, a a /X. As said before, in this type of structures 
both sides of the corrugated metal film are only coupled 
by the term G v that is a function of metal thickness: 
G v = 1/ sin(kW). The diagonal terms e a are related to 
the back and forth bouncing of the EM-fields inside in- 
dentation a: e a — cot(kh a ) for a/0 and eo = cot(fcW) 
for the central slit. Once the values for [E^,E^] are 
calculated, the total transmittance is a function of the 
amplitudes of the E-field at the entrance and exit of the 
central slit: T = Im {E^E^*} /sin(fcVF). 

Also the EM fields in real space can be expressed in 
terms of the set [E^, Eg], For example, in the transmis- 
sion region (z > 0): 



Hy(r) = — Y, E «G(a,r) 



Hoc 



(2) 



all other components of the EM field can be ob- 
tained from H y (r) for the polarization considered. Here 
G(a,r) — ik/2 J <p* a {x)H^\k\xu x — r\)dx. Similar ex- 
pressions can be obtained for the EM-fields in the reflec- 
tion region and within the indentations. 
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III. TRANSMISSION PROPERTIES 

In Ref«i& we analyzed, in a symmetric groove config- 
uration, the conditions for having large transmittance, 
i.e. for having large E-field at the central slit. We 
identified three main resonant mechanisms: slit waveg- 
uide mode excitation (controlled by the metal thickness), 
groove cavity mode (controlled by the groove depths) and 
in-phase groove re-emission processes (that is maximum 
for wavelengths of the order of the period of the array, 
d) . Large enhancements in the transmission for resonant 
wavelengths can be obtained when the three mechanisms 
coincide; this means that if we fix d there is an optimum 
value of the depth of the grooves in which the optical 
transmission is maximum. 

In this section we present results for the total trans- 
mittance T of a system composed of a slit and, in the 
most general case, four periodic finite arrays of grooves. 
We normalize the transmittance to the incident flux in 
the direction normal to the metal surfaces and to the slit 
area. In this way T for point particles would be 1. Even 
with the simplifications introduced, the structure under 
study is still characterized by a large number of geomet- 
rical parameters. In this paper all indentation widths (a) 
are equal and the grooves arrays, when present, have all 
the same number of grooves (N). All results presented 
here are for a = 40nm, N — 10 and metal thickness 
W = 200nm, which are typical values considered previ- 
ously in related experiments (later on we will explain the 
reason for choosing this particular metal thickness) . 

For WDM it is convenient to have at least two wave- 
lengths transmitted primarily through the structure. Let 
us select, arbitrarily and simply for proof of principle pur- 
poses, Ai = 500nm and A2 = 600nm. We consider first 
the even simpler system of a single slit flanked by just 
one groove array in the input side (IL). For the cho- 
sen values of a and N, after an optimization process 
of dh and hi, we find enhanced transmission at Ai for 
di = 470nm,/i£ = 75nm. For this set of parameters, 
groove cavity mode excitation and in-phase groove re- 
emission mechanisms are tuned to appear at Ai. The 
curve T(A) for this case is rendered in Fig. 2 (blue curve), 
showing clearly the transmission resonance at A = 500nm 
due to the presence of grooves (compare this curve with 
the black curve of Fig. 2, showing T(A) for a bare sin- 
gle slit). The green curve of Fig. 2 presents the T(A) for 
a similar system, also with only one groove array next 
to the single slit but this time with lattice parameter 
dji — 560nm and groove depth Hr — 95nm in order 
to have resonant transmission at A = 600nm. Also in 
Fig. 2 (red curve) we show T(A) for the case of a single 
slit flanked by the two groove arrays considered previ- 
ously. This last curve shows that the transmission reso- 
nances associated to each of the arrays survive, although 
somehow modified, the close presence of the other ar- 
ray. The fact that when both arrays are present the 
transmission peaks at A = 500 and A = 600nm have 
the same height is neither casual nor reflects any hidden 
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FIG. 2: Normalized-to-area transmittance plotted as a func- 
tion of the wavelength for a structure denned by a flat exit 
surface and a input surface with grooves at both sides (red 
line), left side (blue line) and right side (green line). Black 
line corresponds to a single slit without corrugations. Inset 
shows the results considering grooves at both sides of the in- 
put surface and different configurations of the output surface. 



property of the system: simply the value of the metal 
thickness (W = 200nm) has been tuned for this to be 
the case. Notice that similar transmission peak heights 
for the selected wavelengths is a bonus for WDM, and 
we just wanted to show that this condition may be ful- 
filled in this kind of structures. The case presented here 
is, however, not necessarily optimal. Probably it is pos- 
sible to make the selected wavelengths to profit from 
slit waveguide resonancesi^iSiiLiS (the third enhancing 
transmission mechanism previously discussed) and get 
even larger transmittances. The inset to Fig. 2 shows 
that in asymmetric structures, as previously found for 
symmetric onesi^, T(A) is mostly independent of the pat- 
terning of the output surface. 



IV. FAR-FIELD ANGULAR DISTRIBUTION OF 
TRANSMITTED LIGHT 

In this section we present results for the far-field angu- 
lar dependence of the radial component of the Poynting 
vector, S r (r), evaluated in the output region. In the 
subwavelength limit, the transmission angular distribu- 
tion of this kind of systems depends only on the out- 
put corrugation, up to a constant that accounts for the 
transmittance at that wavelength 12 . Therefore, it is con- 
venient to define the far-field angular transmission distri- 
bution 1(9) = rS™ r (6) (r > Nd, A), where S nor (X, r) = 
S(X,r)/T(X). Recent studies in the optica l 11 ! 12 and 
microwavei^ regimes have found that 1(9) for a single 
slit flanked by grooves arrays in the output side shows 
beaming effects for certain frequency ranges. In Refii a 
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phenomenological model, based on the hypothesis of the 
existence of a leaky wave running along the surface, was 
enough to explain the beaming phenomena. A subse- 
quent first-principle modeliS supported the basics of the 
phenomenological approach and the origin of the leaky 
surface wave, and was in good agreement with the ex- 
perimental data. Basically, beaming can be understood 
as follows: the EM field leaving the slit couples to a run- 
ning surface wave (which is the result of self-consistent 
wandering back and forth of EM fields at the grooves, 
which act as surface scattering centers) ; this surface wave 
accumulates phase as it travels along the surface and it 
radiates to vacuum at the groove positions. So, for emis- 
sion purposes, the groove array behaves as an effective 
diffraction grating, but illuminated at a (wavelength de- 
pendent) angle, which is related to the phase accumu- 
lated by the surface running wave. Actually, the emis- 
sion intensity is different at the different radiating points 
of this equivalent diffraction grating, and must be cal- 
culated self-consistently from Eq.(l). However the final 
result shows that, at some angles, the extra optical path 
of waves leaving from different grooves cancel the phase 
of the running wave, resulting in a maximum in 1(0). 

Both studiesii*i£ previously referred to considered a 
symmetric distribution of grooves around the central slit 
and, therefore, showed a symmetric I{0). In this paper 
we consider asymmetric structures, showing that beam- 
ing occurs even when a slit is flanked by just a single 
groove array. Also we show that, to a large extend, the 
beam produced by a groove array is not modified if a 
different groove array is placed also at the output metal 
surface but at the other side of the slit. This is illus- 
trated in Fig. 3 where panel (a) shows I(X,6) for a sys- 
tem with a single groove array is placed in the OL con- 
figuration with geometrical parameters — 470nm and 
h,L = 75nm. Notice that, in section 3, we found that 
this set of parameters are optimum in order to have en- 
hanced transmission at Ai = 500nm. Fig. 3a shows that 
this asymmetric structure presents beaming for a range 
of wavelengths and that maximum beaming appears at 
the same resonant wavelength Ai. It is worth noticing 
that even in this asymmetric case, maximum beaming 
occurs at 9 — 0°. A similar representation appears in 
panel (b), this time for an array of grooves in OR con- 
figuration with dji = 560nm and Iir — 95nm which opti- 
mizes both transmission and 1(9 = 0°) for A2 = 600nm. 
Panel (c) renders J(A, 9) when both OL and OR arrays 
previously considered are present in the structure. The 
overall picture, more precisely lines in the (A — 9) plane 
showing large values of /(A, 9) corresponds to the addi- 
tion of beaming lines of panels (a) and (b). However, 
panel (c) also shows that radiation coming from left and 
right groove arrays is not totally independent. Instead, 
there is a small shift in both wavelengths at which max- 
imum beaming appears and the associated angles which 
now may move away from normal direction. This last 
fact will be of paramount importance in order to use this 
kind of structures as wavelength de-multiplexers. 
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FIG. 3: Contour plots of I(\,9) with denned in Fig. 1. 
Three different configurations of the output surface are con- 
sidered, (a) OL with h—75 nm and d—470 nm. (b) OR with 
h—95 nm and d=600 nm. (c) OL+OR with the same geomet- 
rical parameters as in (a) and (b). 
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FIG. 4: Dependence of rS r (0) (see text) on the angle for 
several wavelengths. The result of each wavelength is plotted 
with its corresponding color in the visible spectrum. Notice 
beaming effect present for orange and green light. 



V. WAVELENGTH DE-MULTIPLEXER DEVICE 

The two abilities of these systems described in the 
two previous sections, namely i) the selective transmis- 
sion of two wavelengths when structuring the input side 
of the metallic film (Section 3) and ii) beaming at the 
same two wavelengths in different well-defined directions 
(Section 4) can be combined to propose a new applica- 
tion of single apertures surrounded by finite arrays of 
grooves.: this kind of structures can be used as wave- 
length de- multiplexers of micro- meter dimensions. This 
new functionality is illustrated in Fig. 4 where rS r (0) 
in the far-field region is calculated for different wave- 
lengths for a structure in which four groove arrays are 
present. Geometrical parameters of these arrays are es- 
sentially those considered in Fig. 2 (red curve) and Fig. 3c, 
only that one of the grooves array in the input side has 
been slightly modified in order to obtain enhanced trans- 
mission at A = 530nm. For the two resonant wavelengths 
(A = 530nm and A = 600nm) light emerges in the form 
of very narrow beams presenting very low divergence 
(around 5°) and can be collected at an angle of 6° for 
A = 530nm and 0° for A = 600nm. Note that due to 
the interplay between the left and right grooves arrays in 
the output side, light of wavelength 600nm can also be 
found at a larger collection angle of 14°. For the rest of 
wavelengths (see Fig. 4) optical transmission is extremely 
low and its angular distribution is nearly uniform. 

To further illustrate the resonant nature of this phe- 
nomenon!, in Fig. 5 we show the near-field pictures (E- 
field intensity) in the output region corresponding to the 
transmission processes for the two resonant wavelengths 
(A = 530nm and A = 600nm) and for a non-resonant case 
(A = 700nm). For this wavelength, grooves are not play- 
ing any role; E-field intensity is low and the transmission 
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FIG. 5: Contour plots of the electric held intensity for the pro- 
posed demultiplexing structure. Three different wavelengths 
are considered, (a) A=530 nm. (b)A=600 nm. (c)A=700 nm 
(off-resonance). Red (blue) corresponds to maximum (mini- 
mum) intensity. Black areas show the metallic regions. 



angular distribution is rather uniform. However in the 
first two cases, as clearly seen in panels (a) and (b), light 
is emitted not only by the central slit but also by the 
grooves located at left and right sides of the aperture in 
the output side. As a result of the interference between 
these re-emission processes narrow beams of light emerge 
from the structure and are already present in the near- 
field regime. Fingerprints of the running surface waves 
that are at the origin of the beaming phenomena can also 
be visualized in these pictures. 

In conclusion, we have shown how a single subwave- 
length aperture surrounded by finite arrays of indenta- 
tions may be used as a wavelength de-multiplexer. The 
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two selected wavelengths can be chosen by just modi- 
fying the geometrical parameters defining the structure 
(period of the arrays at left and right of the central slit 
and widths and depths of the indentations). A prac- 
tical advantage of this type of plasmon wavelength de- 
multiplexer device with respect to others based on, for 
example, photonic crystals^ is that light emerging from 
our proposed device forms two collimated beams present- 
ing both a very low divergence. This divergence can also 
be controlled and reduced by increasing the number of 
indentations in the output side. Implementation of this 
device to build up a multi-frequency de- multiplexer can 
be achieved by just adding up other single apertures and 
the corresponding indentations with the appropriate ge- 
ometrical parameters to the structure analyzed in this 



communication. The fact that input and output corru- 
gations seem to act independently and also that the in- 
terplay between left and right indentations is weak allows 
us to state that this type of structured metals could be 
used as building-blocks of several opto-electronic devices 
of reduced dimensions. 
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